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Abstract 

In this paper we investigate the conformation statistics of a Gaussian chain 
embedded in a medium of finite size, in the presence of quenched random ob- 
stacles. The similarities and differences between the case of random obstacles 
and the case of a Gaussian random potential are elucidated. The connection 
with the density of states of electrons in a metal with random repulsive im- 
purities of finite range is discussed. We also interpret the results obtained in 

some previous numerical simulations. 

PACS numbers: 36.20.Ey, 05.40-a, 75.10.Nr, 64.60.Cn 



1 



Typeset using REVTgX 



I. INTRODUCTION 



The behavior of polymer chains in random media is a well studied problem both theoret- 
ically and experimentally [0-12] and has applications in diverse fields. For polymers, 



the interest arises when the chains are confined inside an intertwined gel network [12], and 
perhaps inside porous materials and membranes P-P"Hj. Furthermore, the problem is related 
to the statistical mechanics of a quantum particle in a random potential ||13|| , the behavior 
of flux lines in superconductors in the presence of columnar defects [13,13], and the problem 
of diffusion in a random catalytic environment 

In this paper we will study the static properties of a Gaussian polymer chain, without 
excluded volume interactions, that is confined in a quenched random medium. Experimen- 
tally, a polymer in a specific solvent is known to obey Gaussian statistics at the so called 
temperature-where the long range self-avoiding interactions are effectively screened. The 
term quenched refers to the fact that the random medium is frozen and thus does not ther- 
malize with the active degrees of freedom-in this case the polymer chain. We are interested 
in the properties of the polymer, such as the free energy and radius of gyration (or alter- 
natively the end-to-end distance), that are averaged with the appropriate Wiener measure 
over all possible chain conformations in a given realization of the random medium, with a 
final average taken over all possible configurations of the random medium. The nature of 
the random environment is crucial in this problem, and so it is important to distinguish the 
following two important cases that have been discussed in the literature: 

1. A Gaussian random potential with short range correlations. 

2. Random obstacles which prevent the chain from visiting certain sites. 

Numerical simulations performed in three dimensions were restricted, to our knowledge, 
only to the case of random obstacles |1]|||6[]. Also, a numerical investigation based on the 
mapping to the Schrodinger equation in one dimension was performed for the Gaussian 
potential case ||. On the other hand, extensive analytical work using the replica 
variational approach and Flory type free energy arguments, has been done for the case of a 
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Gaussian random potential . The case of a bounded (saturated) random potential 

was also addressed in ||. It was not clear to us to what extent these theoretical 
investigations could be applied to the case of infinitely strong random obstacles placed 
randomly in the medium, as simulated numerically. 

In this paper we investigate analytically, for the first time, the case of infinitely strong, 
randomly placed obstacles case. We point out similarities and differences with the case of a 
Gaussian random potential, and also the case of a saturated potential. We will assume that 
the obstacles are infinitely strong-they totally exclude the chain from visiting a given site 
occupied by an obstacle. Each obstacle is taken to be a block of volume a d , where d is the 
number of spatial dimensions and where a is the linear dimension of the block. We take for 
simplicity the polymer bond length b to be approximately equal to a. Thus, a will be the 
small length scale in the problem, and we will measure all distances in units of a. However, 
in the next section we will sometimes keep a explicitly in order to omit terms of higher order 
of smallness. The obstacles are placed on the sites of a cubic lattice with lattice spacing a. 
We denote by x the probability that any given lattice site is occupied by an obstacle (block). 
Our main results will concern the case of small x, in particular x < x c , where x c refers to the 
percolation threshold (x c = 0.3116 for a cubic lattice in d = 3), but we will also comment 
on the case of a larger concentration of obstacles. We denote by V the total volume of the 
system. 

For an uncorrelated Gaussian random potential, it was argued using qualitative argu- 
ments in Refs. |3|,[|] that a very long Gaussian chain will typically curl up in some small 
region of low average potential. The polymer chain is said to be localized, and for long 
chains the end-to-end distance (Rf) becomes independent of chain length and scales like 

R F (x{g\nV)- l/{A - d \ (1) 

with g being the strength of the disorder (the random potential satisfies (U(x)U(x')) = 
gS(x — x')). The depth of the well entrapping the chain is approximately 

U min ~ -(glnV) 2 ^. (2) 
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These results were also obtained by the replica method in Ref. [0], and rederived using a 
mapping to a quantum particle's localization in Ref. ||. For very short chains the end-to- 
end distance scales diffusively (R F ~ L), and it saturates at the Rf value quoted above 
for large L. Notice, that in the infinite volume limit, the chain completely collapses. This 
results from the fact that the depth of the potential is unbounded from below, and the chain 
is always able to find with reasonable probability a deep enough narrow potential well to 
occupy, overcoming its tendency to swell due to the entropy of confinement. The collapse 
of the chain in the infinite volume limit agrees with the results for a chain in an annealed 
potential, since the ability of a chain to scan all of space for a favorable environment is 
equivalent to the random potential adapting itself to the chain configuration. 

To review briefly the argument leading to Eqs. (|], 0) using localization theory || we 
recall that the density of states for a particle in an uncorrelated Gaussian random potential 
is given by |L6] 



p(E) = j^ex V (-B\E\ s ), (3) 

with 5 = (4 — d)/2 and B oc 1/g, where g is the strength of the disorder. This result is valid 
for an infinite volume. In a finite volume V the energy will be bounded from below. We can 
estimate the lowest energy E c from the tail of the distribution: 

dEp(E) ~ 1/V, (4) 



In 1/5 



which leads to 

The width of the ground state wave function (localization length) is given by 

4~|^r 1/2 ~(^lnV)- 1/(25) . (6) 

The mapping from a quantum particle of mass m, at a finite temperature 1/(3, to a polymer 
chain, is given by 

% -> T, h(3^L, m ^ dT/b 2 . (7) 
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It can be shown that the ground state width £ c is proportional to the end-to-end distance 
Rp of a chain that is situated in a deep minimum in the volume V, whose depth is given 
by E c ~ U m i n . Using the above mapping the density matrix for a quantum particle at finite 
temperature |T7|] corresponds to the partition sum (Green's function) of a Gaussian polymer 
chain [IS|. 



II. RANDOM OBSTACLES 

For the case of infinitely strong randomly placed obstacles, the potential energy of the 
polymer chain is always zero. Hence, the free energy of the polymer will be F = —TS, since 
E — 0, and the statistics of the polymer will be dictated only by entropic effects. As the 
volume of the system tends to infinity there is a chance to find very large lacunae free of 
obstacles. Thus, in the limit of large volume we do not expect the polymer to collapse, but 
rather to inflate with increasing number of monomers (L). We will now analyze the behavior 
of a polymer in an environment consisting of random obstacles and find that there are three 
different phases as a function of the volume of the system. In the subsequent analysis we 
will always assume that L is very large. 

In order to estimate the average chain properties we first assume that a very long polymer 
chain will attain an approximately spherical shape of radius R (we will discuss this spherical 
droplet approximation in the sequel). Now, let us coarse grain the volume V into subregions 
of volume v ~ R d , and assume that the polymer is confined to one of these regions. Each of 
these coarse grained subregions will contain a different number of obstacles, and the chain 
will reside in that region with the lowest number of obstacles which can be found in the 
finite volume V in order to minimize its free energy. We will estimate R by writing the free 
energy of the polymer as a function of both R and the coarse grained volume fraction of 
obstacles in a given subregion (which also depends on R), and minimizing it accordingly. 
First, let us assume that there are no obstacles present inside this spherical region. The 
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entropy of a chain confined in this "cavity" is of the form 

S = L\n(z) - a^, 

where z is the number of nearest neighbors, e.g. 2d for a cubic lattice in d dimensions, and 



the second term is the entropy of confinement ||19|| . Here, a is a numerical constant. The 
free energy will then be F = —TS. In the following we choose T = 1 for simplicity since 
the temperature does not play any significant role with respect to the results. 

In order to proceed and estimate the entropy change due to obstacles inside the volume v, 
we use the mapping from the polymer to a quantum particle as discussed in the introduction. 
The free energy per monomer of the polymer when L is very large corresponds to the 
ground state energy of a quantum particle in a cavity of radius R. This is known, in three 
dimensions, to be equal to E = (h 2 /2m)7r 2 / R 2 , in agreement with the expression above (up 
to an unimportant additive constant). In d dimensions the energy is still proportional to 
1/R 2 but the prefactor is different. Suppose that there is a spherical obstacle of radius a 
inside the sphere. If the obstacle is at the center of the sphere the Schrodinger equation is 
exactly solvable and the ground state in d = 3 is given by 

7T (r-.fi) 



sm 



*( r ) = C "~~ R ~ a , a<r <R (8) 



and \l/(r) =0 otherwise. This will correspond to an energy of 

m it 2 iT 2 a 

¥ Eq = 2R 2 + 1v + --- (9) 

where the corrections vanish faster than a as a — > 0. If on the other hand the obstacle is not 
in the center of the sphere we could only find an approximate solution of the Schrodinger 
equation which can be used to give an upper bound to the ground state, which may be exact 
to leading order in a (see Appendix). The ground state energy becomes 

m 7r 2 7r 2 a / R . tcRq^ 2 



J? E °=2IP + 1P {^ Sm lf) + — (10) 
where Rq is the distance of the center of the obstacle from the center of the sphere. One 
can see that the factor in parenthesis approaches 1 as Rq — > 0, and vanishes as Rq — > R. 



Notice, that for the analysis above we have treated the obstacles as spherical in shape as 
opposed to a square block. However, this difference should only amount to an unimportant 
numerical factor. 

Using the mapping from the quantum particle to the polymer as given by Eq. ([/]), we 
find 

(«) 

where F is the free energy of the chain. We have used the fact that T = 1 and b = a for 
the bond length. (This correspondence is up to an additive constant that is given by the 
entropy of a chain in free space). Now, if the position of the obstacle is random we expect 
that we should average the energy over all the possible locations of the obstacle within the 
volume of the sphere to obtain 

3F _ vr 2 3a 

~ + ^ + ••• • ( 12 ) 

Let us denote by x the volume fraction occupied by the random obstacles within the spherical 
volume v. Then the number of obstacles inside this volume will be ^-(R/a) 3 x. If x is small 
the energy due to several obstacles will be approximately equal to the sum of the individual 
energies. The deviation from this rule becomes important only if the obstacles are very 
close or touching each other, and for small x the number of such configurations is very small 
compared with configurations where the obstacles don't touch. In any case interactions 
among obstacles are at least of 0(x 2 ). Thus, the free energy of a polymer chain confined to 
a volume v, with a volume fraction x occupied by obstacles, is 

F ^ 27r ~ 

L = 6(W + ^ + "-- (13) 
The important conclusion is that the term proportional to x is independent of R. The nu- 
merical prefactors are really of no importance to us. The small x approximation is enhanced 
by the fact that not only do we assume that x is small (x < x c ), but the chain will be found 
to settle in regions where the obstacle concentration is smaller than average. In two dimen- 
sions we find that the first term in Eq. (|T2| ) is proportional to 1/R 2 and the second term to 



1/(R 2 \ \n(a/R)\) (see Appendix). In one dimension the situation is totally different since an 
obstacle will divide the volume v into two disjoint regions. Thus, all our conclusions apply 
only above two dimensions, where the second term in Eq. ([12]) is proportional to a d ~ 2 R~ d . 
When multiplying by the number of obstacles ~ (R/a) d x one gets an /^-independent result 
for the second term in Eq. (0). Two dimensions is a borderline case where some subtleties 
may arise. In the following we will measure all distances in units of a so we put a = 1. 

We now proceed to study the chain statistics by considering the fluctuations of the 
volume fraction x. If x is the probability that a lattice site is occupied by an obstacle, then 
the coarse grained volume fraction x is distributed according to the binomial probability 
distribution b(vx;v,x). The notation 

b(k;n,p)=(fjp k q n - k , (14) 

stands for the probability that n Bernoulli trials with probabilities p for success and q = 1 — p 
for failure result in k successes and n — k failures ||20|| . For x — it gives the so called Lifshits 
probability to find an empty region of volume v free of obstacles. With this probability there 
will be associated an "entropy" which will be its logarithm. 

Using the results above we can start to discuss the statistics of a polymer in an infinite 
volume (V — > oo)— the so called annealed result. Assuming that the chain takes on a roughly 
spherical configuration of volume v ~ R d , the free energy will read 

F(R, x) = -L \n(z) + -^ + Lx- \n[b(R d x; R d , x)] . (15) 

This free energy has to be minimized both with respect to x, and to R, since the chain is free 
to move and find the most favorable values for these parameters. The most favorable value 
of x for large L and for an infinite volume is x = 0, since x is not allowed to be negative. 
Using the fact that 

b(Q;v,x) = (l-x) v , (16) 
we find that the expression for the free energy becomes 

F(R) = -L\n(z) + -^-R d \n(l-x). (17) 



This free energy has now to be minimized with respect to R to yield 

f L \ V( d + 2 ) 

Rm.annealed ~ I 7~J 7Z 77 I • (^) 

V|ln(l 

Thus the size of the chain grows with L, but with an exponent smaller than 1/2, the free 
chain exponent. 

So far we discussed the case of an infinite volume V. In a finite volume we find that the 
so called quenched and annealed case differ, at least when the volume is not too big. We 
actually find that there are three regions as a function of the size of the system volume V. 
First, if V < Vi — exp(x~( d ~ 2 ^ 2 / '(1 — x)), it is unlikely for a chain of volume v ~ R d to find 
a region which is totally free of obstacles. Thus x does not vanish in this regime. To proceed 
further we must use an approximation to the binomial distribution b(vx; v,x). 

If v is not too small we can approximate the binomial distribution by a normal distribu- 
tion gg 

b(vx; v, x) « (2ttvx(1 - x))~ l/2 exp f _ v ( x ~ x ) \ _ / ig \ 

\ 2x(l — x) J 

This approximation is good provided vx ^> 1 and v(l — x) ^> 1. We will verify below that 
these conditions are indeed met in our case when x is small. 

In a finite volume V, the lowest expected value of x, to be denoted by x m , can be found 
from the tail of the distribution 



o 



I -^f) - ?• (20) 



which gives 



xy InV , , 



The free energy becomes 



Fj(R) = -L ln(z) + ±- 2 +Lx-L\l (22) 

The last term in Eq. (|T5| ) is missing since it is negligible for large L when R is independent 
of L. Minimizing F(R) with respect to R we find 

R mI ~ {xy\nV)- 1/{A - d) (23) 
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and 

x mI = x - {xy\nVf /{i ~ d) , (24) 

where we put y = 1 — x. The result for the radius of gyration of the chain, as represented 
by R m i is the same result as for the case of the Gaussian distributed random potential, but 
with the strength g replaced by x(l — x). The polymer in this case is localized and its size 
is independent of L for large L. 

As V grows R m decreases until eventually x m vanishes. This happens when V = Vi — 
exp(x~( d ~ 2 ^ 2 y~ 1 ) . For V > Vi, R m is no longer given by Eq. (0), but rather by the solution 
of x m u = 0. It is the largest region free of obstacles expected to be found in a volume V. 
Rather than using the normal approximation we can estimate R m directly from the relation 

(1 - x) v ~ v/V, (25) 

with v ~ R^. Solving for R m we obtain 

R^, (26) 



ln(l — x) 

The polymer is still localized but the dependence on x and on InV has changed. In this 
region which we call region II the free energy is given by 

(1 V N — 
• (27) 

where some undetermined constant a is introduced for later convenience. 

As V grows in region II, R m ii continues to grow until it reaches the annealed value given 
in Eq. (JTH). This happens when 

V = V 2 ~ exp ( x V(d + 2) L d/ {d+ 2)^ (2g) 

to leading order in x, which is enormous for large L. For V > V2 we have the third region 
in which R m in = R m ,anneaied and it grows like L 1/(a!+2) . 

Since in region I we have used the normal approximation to the binomial distribution, 
we should check a posteriori if the condition vx ^> 1 is met. Since V < V\ we find 

C d R d mI x » C d x-^ > C d x-^l\ (29) 
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Here C d = 2-n d l 2 / (dT(d/2)) is the volume of a ci-dimensional sphere of unit radius. In 
d = 3, for example, C3 = Att/3. Thus we see that for d > 2, vx ^> 1 provided a; < 1. 
The minimum value is attained for V = Vi, and is significantly larger when V < Vi. It is 
interesting to notice that in d = 3 for example, even if x = 0.3, we find that C3/ ' \fx = 7.6 
which is considered large enough for the validity of the normal approximation ( a value of 5 
is usually considered sufficient). 

Actually, the normal approximation to the binomial distribution b(k; n,p) is not entirely 
acccurate even if the conditions np ^> 1 and nq 3> 1 are met, if k is far from the center, i.e. 
if \k — np\ 3 /n 2 > 1 ||20|. In our case we need that v\x m — x\ 3 < 1 to be satisfied. In region I 



we have 

C d R d mI \x mI - x| 3 < C d (xy\n V^'^ = C d x^l\ (30) 

where we have used Eq. ([24]) and the value of V\. For d = 3 this is smaller than one even 
for x = 0.3. 

The behavior in region II can also be deduced from known results of the density of states 
for a quantum particle in the presence of obstacles (repulsive impurities). In that case [TH 
the density of states is given by (when the obstacles are placed on a lattice) 

p(E) ~ exp(-c| ln(l - x)\E- d/2 ), E > (31) 

with c being some dimension dependent constant and x is the density of impurities. Note 
that p(E) vanishes for E < 0. We can estimate the lowest energy in a finite volume V from 
the integral 

dEp(E) ~ 1/V, (32) 







and find 



InV ^ 2/,/ 



E *~ii4rw ] • (33) 



and thus the localization length is given by 



InV ^ 
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which coincides with R m ii above. 

To conclude this section we make some remarks on the validity of the spherical droplet 
approximation. The shape of a long polymer chain is determined by the regions of the 
random medium that have a lower than average number of obstacles. For V > V\ these 
regions are essentially free of obstacles. The probability of finding such empty regions 
depends only on its volume and not its shape. However given regions of varying shapes 
and equal volumes, it will be entropically more favorable for a long polymer chain to reside 
in a region whose shape is closest to a sphere. This is because the confinement entropy is 
maximized for a sphere over other shapes of the same volume. The argument is equivalent to 



that proposed by Lifshits |T6| in the context of electron localization and is shown rigorously 
by Luttinger ||21|| . For V > Vi the relevant regions contain a small number of obstacles 
but we believe that the same argument should roughly hold and deviations from a spherical 
shape will be small or irrelevant. 



III. COMPARISON WITH NUMERICAL SIMULATIONS 

We have compared our results from the last section with numerical simulations performed 
by Dayantis et al. ||, and also comment on the relation to earlier simulations done by 
Baumgartner and Muthukumar [p]]. Dayantis et al. carried out simulations of free chains 
(random-flight walks) confined to cubes of various linear dimensions 6 — 20, in units of the 
lattice constant. These chains can intersect freely and lie on a cubic lattice. They introduced 
random obstacles with concentrations x = 0, 0.1, 0.2 and 0.3. The length of the chains vary 
between 18 — 98 steps. They also simulated self- avoiding chains that we will not discuss 
here. They measured the quenched entropy, the end-to-end distance, and also the radius 
of gyration which is a closely related quantity. Unfortunately, these authors did not have 
a theoretical framework to analyze their data, and thus could not make it collapse in any 
meaningful way. We show below how it is possible to fit the data nicely to our analytical 
results. 

Even for x = 0.1, the value that we get for Vi is about 33 which is an order of magnitude 
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FIG. 1. A plot of ln(— S/L + In 6) vs. ln(ln V/| ln(l — x)\). The labels are marked according to 
the chain length. 

smaller than the the smallest volume used in their simulation, which is 216 for a cube of side 
6. Hence we expect to be in region II. To check the agreement with our analytical results 
we show in Figure 1 a plot of ln(— S/ L + In 6) vs. ln(ln V/| ln(l — x) |) where 5* is the entropy 
measured in the simulations and V = B 3 for a box of side B. Recall that F = —S and 
Eq. (p7|) predicts a straight line with slope —2/3. The best fit is obtained for a slope of 
—0.72 ± 0.05, which is in excellent agreement with our analytical results in region II. 

In order to analyze the simulation results for the end-to-end distance and radius of 
gyration we have to introduce some additional compensation for the results obtained in the 
previous section. First we must realize that Eq. (|26|) is valid only when the number of steps 
(monomers) is very large. In the simulations they used chains of varying lengths whose size 
did not yet reach asymptotia. Hence, we introduce a correction factor 

R m (L) = R m {l - expi-L/Rl)) 1 / 2 = R m h{^L/R m ) , (35) 

which interpolates between the size of a free chain as L — > and the value of R m from 
Eq. ( |26"D as L — > oo. 

The second correction we have to implement arises when the expected value of the chain 
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is not much smaller than the size of the confining box. Even for a free chain confined to 
a box of side B with no obstacles present, the end-to-end distance is not simply R = L 1//2 
for L 1 / 2 < B and R = B for larger L. We have to take into account the fact that the 
length of the chain has a Gaussian distribution about its expected value, and the tail of the 
Gaussian is cut off by the presence of the box (this is for the absorbing boundary conditions 
that is used in the simulations). Thus, for the case of no obstacles (x = 0), The measured 
end-to-end distance should approximately be 

Rl = J* dR&exp{~)/ J B dR exp(-^), (36) 

which gives 

R c = y/Lf 2 {B/VZ) with 

/2W ^( i -/!sR7i)^ 2/2) ) I/2 - (37) 

This indeed gives good agreement with the measured values in the no obstacle case. For the 
obstacle case we thus have to introduce these two corrections in subsequent order: 

R m,corrected R m h{^fl/R m )h - j—— . . , (38) 

\RmJlW L / R rn> / 

where R m = R m ii as given by Eq. fl2"6"|). (A constant of proportionality of 1.8 has been 
introduced on the rhs of Eq. fl26|) to obtain a good fit). In Figures 2 and 3 we show 
a comparison of the simulation results for the end-to-end distance and for the radius of 
gyration with the calculated results as given by Eq. ( p6|) and Eq. (|38|). The agreement 
seems remarkable, especially for the end-to-end distance, where all the data collapses to a 
straight line with a slope close to 1. 

Dayantis et al. emphasize that they did not consider concentrations of obstacles above 
the percolation threshold, which is at x c = 0.3116 for a simple cubic lattice. The reason 
is that above the percolation threshold the medium of random obstacles begins to form 
disconnected islands free of obstacles. Thus, in their simulation the polymer chain will only 
sample a limited fraction of the volume available. What happens is that effectively the 
volume available for the chain is not the total volume of the cube but rather the volume 
of the disconnected region it occupies. For most realizations of the random medium this 

14 




calculated 



FIG. 2. A plot of the observed vs. calculated end-to-end distance 




FIG. 3. A plot of the observed vs. calculated radius of gyration 
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effective volume will be smaller than the value Vi, which is the limit of region I of the 
last section. In that case one expects the end-to-end distance to scale like x _1 as given 
in Eq. fl23|) instead of like a; -1 / 3 as given by Eq. fl26|). Baumgartner and Muthukumar's 
simulation was for both below the percolation threshold and also above it (x = 0.4 and 0.5). 
However, they only estimate the exponent above the percolation threshold, and find it to be 
about —1. They do not estimate the exponent for x below the percolation threshold, which 
appears from their data to scale with a much smaller exponent. Thus, it seems likely that 
the reason these authors report a behavior corresponding to region I, even though their box 
is quite large, is because the effective volume is small for the cases for which they exceed 
the percolation threshold. 



IV. SUMMARY AND DISCUSSION 

In this paper we have considered the effect of random obstacles on the behavior of a 
free Gaussian chain. We have seen that in the presence of infinitely strong obstacles, that 
exclude the chain from visiting randomly chosen sites, there are three possible behaviors of 
the end-to-end distance as a function of L (the number of monomers). The three possible 
regimes depend on the total volume of the system. If the volume of the system is smaller 
than Vi — exp(x^ ( - d ^ 2 ^ 2 ), where x is the probability that a lattice site is occupied, then the 
polymer is localized, as in the case of a Gaussian random potential, and Ri ~ (x In V) 1 '^~ d ' . 
For Vi < V < V2, where V2 ~ exp (x 2 '( d+2 ^L d '( d+2 n, the polymer size is given by Rn ~ 
[xj In V) _1//d . Finally, for V > V2, the polymer behaves the same way as for an annealed 
potential, i.e. Rui ~ (L/x) l ^ d+2 \ We have displayed the results to leading order in x for 
small x. These results are valid only when the average volume fraction of the obstacles (x) 
is smaller than the percolation threshold. When x is bigger than the percolation threshold 
we expect that the system breaks into independent domains whose volume is independent 
of, and generally much smaller than, the total volume of the system. 

We were able to fit nicely the simulations of Dayantis et al. U to the analytic behavior in 
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region II. The results of Baumgartner and Muthukumar for x > x c , presumably correspond 
to region I. Region III is very difficult to be seen in simulations since V2 is huge for large 
L. However such behavior, which coincides with the annealed case have been obtained by 
Chandler et al. || in their simulations of annealed random obstacles. 

The behavior in the three regions has some similarities to the case of a saturated random 
potential as discussed by Cates and Ball [|J. However, there are also significant differences. 
The saturated potential discussed by these authors concerned the case where the random 
potential can take two different values iy 7 ^ with probability 1/2 each. They considered 
the case when g is small. In that case, when the chain occupies an approximately spherical 
volume R d , its energy is estimated by the average potential in that region times the number 
of monomers L of the chain. This leads to a coarse grained renormalized potential with 
a Gaussian probability distribution that is valid provided the predicted size of the chain 
is much larger than the largest cavity that is free of obstacles. On the other hand, in the 
case of infinitely strong obstacles the free energy of a chain occupying a volume R d has a 
completely entropic origin. The energy of the chain remains zero as it is completely excluded 
from the regions occupied by obstacles. In the case of a saturated potential the behavior 
in region I is given by R$i ~ (gin V)" 1 ^ 4- ^, so in this region g plays a similar role to x 
(the subscript S refers to the saturated potential case) . However, the volume of the system 
beyond which this behavior is no longer valid is given by Vsi — exp(g~ d / 4 ), which is quite 
different from its value in the strong obstacle case. Here d = 2 is no longer a lower critical 
dimension, and the result, which has a very different <i-dependence, is valid down to (and 
including) one dimension. For Vsi< V < V52, one has Rsn ~ (lnV) 1//d , independent of g. 
Also Vs2 ~ exp (L d/ '( d+2 )) , independent of g. Above V52 the annealed result Rsin ~ L 1 /^ 2 ) 
applies. 

As compared to the unsaturated gaussian random potential, the two models coincide 
only in region I when V < V±. For V — > 00 the polymer collapses to a point in the Gaussian 
potential case as it can find a very deep and narrow potential well, whereas in the random 
obstacle case the polymer swells as L l ^ d+2 ^ with growing L as it can find large lacunae free 
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of obstacles. Our analysis of the simulations of Dayantis et al. || shows clearly that their 
results do not conform to the behavior dictated by the Gaussian model but are described 
well by the expected behavior in the intermediate region II. 
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APPENDIX: 

Consider a spherical cavity of radius R about the origin. Inside the cavity there is one 
obstacle of radius a, centered at location Ro. The potential within the obstacle is infinite 
and thus the wave function has to vanish in that region. We will assume that a <C R. Let 
us define the unperturbed wave function 

/ \ A 7r|r| . , . 

^o(r) = ^sin-^. (Al) 

This is the ground state solution of the Schrodinger equation in the absence of the obstacle. 
Consider now the following trial wave function: 



with 



0, < |r - R | < a 

*o(r)*i(r-Ro), a < |r - R | < r m (A2) 

*o(r), r m <|r-R | && Irl < R 



, m 1 . / 7r(|r| - a) , 

ir(r m -a) \ W \ R 



sm, R 



Here r m is determined by the condition 

V* 1 (r)| H=rro = 0, (A4) 



18 



which leads to the condition 



tan " (r ^ a) =^. (A5) 



In the limit of small a this leads to 

'4) i? 2 /V/ 3 . (A6) 



' m — 



3X l/3 



7l 2 



Thus the volume of the region around the obstacle in which the wave function deviates from 
\I>o is about 4R 2 a/ir, which vanishes like a as a — > , not like a 3 . In the following we 
will refer to this volume as V a . The choice of r m is designed to insure that not only the 
wave function but also its gradient are continuous across the seam |r — R | = r m . Thus As 
|r — R | — > r m from below 

Vtf = ^iW + ^oWi -> W , (A7) 

since Wi -> and *i -> 1. 

Using this wave function an upper bound on the energy is given by 

" JW)^ • (A8) 

Evaluating the rhs we find 

" + /*§(r)dV + 0( ° } ' (A9) 

where in the numerator we only integrate over the volume V a about the obstacle. Now, 
inside V a 

V 2 ^(r) = *i(r - R )V 2 ^ (r) + ^ (r)V 2 ^!(r - R ) + W (r) • Wi(r - R ). (A10) 
We argue that 

J W (r)-V$i(r-R ), 

vanishes faster than a in the limit a — > 0, because of the angular integration: V^oir) points 
in a fixed direction away from the origin at 0, whereas V\l/i(r — R ) points in a varying 
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radial direction away from the center at Ro, and thus the angular integration will yield 
cos 9 dcosO = 0. We have actually verified that the contribution of this term is 0(a 5//3 ), 
hence negligible. Since -(1/2) V 2 \t> (r) = (7r 2 /2i? 2 )\l/o(r), the first term in Eq. ([MUD exactly 
cancels the second term of the integral in the numerator of Eq. (|A9). Thus we are left with 



m p vr 2 / gg(r)*i(r ~ R )(-|V 2 )^i(r - R )dV a 

W E °-2& + jn(r)dV • (AU) 

At this point we can approximate ^o( r ) by its value at r = R , and take out of the integral. 
Also, 

(-iv 2 )^ - Ro) = (Tr 2 ^ 2 )^ - Ro). 

Thus 



m vr 2 7r 2 . [ ty?(r)dV a 

/•% + »*»7ipF' (A12) 

Evaluating the integrals and using the estimate for r m as given by Eq. (|A6|) , we arrive at 
the final result 

777 7T 2 7T 2 a / R TtRq\ 2 



ft 2 ~~ W? ~R? \kR~q ~ R - 

as given in section II, where Rq is the magnitude of Ro- Although we have used an approx- 
imate wave function that can give only an upper bound on the correction, it appears likely 
that the correction might be exact to leading order in a. Support for this comes from the 
fact that it coincides with the exact answer in the limit Rq = 0. 

We now consider briefly the two dimensional case. In this case we consider only the case 
where the obstacle is in the middle of the circular region. Without the obstacle present, the 
ground state solution to the Schrodinger equation is given by 

2.405r\ 



*o(r) = AIo ( — — ) , (A14) 



where Jo is the spherical Bessel function and xq ~2.405 is its first zero. This leads to 
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With an obstacle of radius r present at the center, we look for a wave function of the form 

{AJ (At) + BN (Ax) , r> a 

(A16) 
0, r < a 

where N is Neumann's function, sometimes denoted by Y . The variable k is determined 
by the requirement that ^o(r) = at r = a and r = R. This leads to the condition 

Jo (kR) J (ka) 



N (kR) N (ka) 



(A17) 



Thus 



3 (kR) l + 0(k 2 a 2 " 



and we find 



N (kR) (2/7r)(ln(fca) + 7 - In 2) + 0(k 2 a 2 ) ' 



(A18) 



* « + "rjimTT' (M9) 



i? #|m|| 



with 



TrNoOro) , Aon x 

" = 2 JiM ' (A20) 

Here J is the derivative of J with respect to its argument. It thus follows that 

W E ° ~ 2R~ 2 + Wfh[f\- (A21) 
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